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SUMMARY 

The  energy  transfer  from  the  uniform  flow  outside  a  boundary 
layer  to  a  transverse  surface  wave  at  the  boundary  is  calculated  on 
the  hypothesis  that  the  boundary  layer  may  be  represented  by  an 
inviscid,  approximately  parallel  shear  flow.  This  energy  transfer 
is  found  to  consist  of  two  components;  the  first  is  similar  to  that 
found  previously  for  supersonic  panel  flutter  in  the  absence  of  a 
boundary  layer  and  is  relatively  diminished  by  the  presence  of  the 
boundary  layer;  the  second  is  intrinsic  to  the  shear  flow  and  iB 
present  whenever  the  surface  wave  speed  is  smaller  than  the  free 
stream  speed  (whether  subsonic  or  supersonic).  Approximate 
solutions  to  the  differential  equations  for  the^dlsturbed  motion  of 
the  boundary  layer  are  establtehed  and  compared  with  more  exact 
jauiweeieal  solutions  ■  ^The  application  of  the  results  to  actual  flutter 
calculations  is  discussed  and  an  example  of  a  pressurized  cylinder 
considered.  It  is  concluded  that  the  presence  of  a  typical  boundary 
layer  may  reduce  the  degree  of  instability  for  supersonic  flutter  of 
a  long,  monocoque  shell  by  an  order  of  magnitude,  thereby 
providing  a  possible  explanation  of  discrepancies  between  earlier 
theories  (in  which  boundary  layer  effects  were  neglected)  and 
observation.  (  j 
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1.  INTRODUCTION 

The  dynamic  instability  (flutter)  of  an  infinite  panel  exposed  to  a 
supersonic  air  stream  has  been  discussed  in  references  1  and  2.  It  was 
found  there  that  the  flutter  speed  was  approximately  +  (Vo)  ,  where 
denotes  sonic  speed  and  VQ  the  speed  of  a  transverse  wave  on  the  free 
panel,  and  that  the  wavelengths  for  monocoque  structures  were  sufficiently 
small  to  justify  the  neglect  of  panel  curvature  on  the  aerodynamic  (but  not 
the  structural)  forces  and  to  justify  the  assumption  of  an  infinite  panel.  The 
smallness  of  the  wavelength  implies  that  the  effects  of  a  relatively  thick 
boundary  layer  may  be  important,  and  the  following  analysis  is  intended  to 
establish  the  qualitative  nature  of  boundary  layer  effects  and  to  provide 
approximations  to  their  magnitude. 

We  adopt  as  a  model  for  our  calculations  an  inviscid,  parallel  shear 
flow  over  an  infinite,  approximately  plane  panel.  Assuming  a  periodic 
surface  wave  displacement  of  the  panel,  we  will  establish  the  equations  of 
motion  for  this  shear  flow  and  the  resulting  pressure  on  the  boundary  in 
Section  2,  after  which  the  eigenvalue  problem  for  the  panel  will  be  established 
in  Section  3,  We  then  will  consider,  in  Section  4,  the  general  nature  of  the 
possible  instability  and  show  that,  in  addition  to  the  essentially  supersonic 
instability  of  reference  2,  there  also  appears  a  new  type  of  instability 
associated  with  energy  transfer  from  the  shear  flow;  this  latter  mechanism 
is  closely  related  to  that  which  enters  the  boundary  layer  stability  problem. 

The  boundary  value  problem  considered  in  Section  2  leads  to  a  second- 
order  differential  equation  that  has  a  regular  singularity  in  the  range  of 
integration  if  the  surface  wave  speed  is  less  than  the  air  speed  at  the  outer 
(free  stream)  edge  of  the  boundary  layer.  We  consider,  in  Section  5,  the 
nature  of  solutions  to  this  differential  equation  in  the  neighborhood  of  its 
singularity,  its  solution  by  expansion  in  powers  of  the  ratio  of  boundary 
layer  thickness  to  wavelength,  and  its  solution  by  numerical  integration. 

The  application  of  these  results  to  actual  flutter  analyses  is  discussed  in 
Sections  6  and  7. 
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We  emphasise  that  the  model  on  which  the  present  analysis  is  based 
is  highly  idealised,  and  experimental  confirmation  of  the  results  is  most 
desirable;  moreover,  even  though  the  predicted  conditions  for  the  initiation 
of  flutter  should  prove  correct,  nonlinear  effects  must  be  decisive  in 
determining  the  severity  of  such  flutter.  In  the  absence  of  experimental 
confirmation  and  at  least  a  qualitative  understanding  of  these  nonlinear 
effects  (especially  with  reference  to  structural  fatigue),  our  results  must 
be  accepted  with  considerable  reserve  and  should  not  be  regarded  as  adequate 
for  major  design  decisions.  It  appears,  nevertheless,  that  typical  boundary 
layers  may  reduce  the  negative  damping  factor- -i.  e.  ,  the  degree  of 
instability- -by  an  order  of  magnitude  (this  statement  applies  only  to  long, 
monocoque  cylinders,  for  which  the  critical  wavelength  and  boundary  layer 
thickness  are  of  the  same  order  of  magnitude),  thereby  providing  a  possible 
explanation  of  previous  discrepancies  between  theory  and  observation. 
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2.  EQUATIONS  OF  MOTION- BOUNDARY  LAYER 


We  require  the  equation#  of  motion  that  govern  the  disturbed  flow  produced 
in  a  prescribed  boundary  layer  by  the  surface  wave  displacement**1 


C  =  ;  (x-  Vt,y)  =  exp 


{ - 1 


(x  -  Vt)  cos  a  +  y  sin  a 


] 


(2.1) 


of  the  boundary  from  its  equilibrium  position  z  =  zq;  x  and  y,  together  with  z, 
are  Cartesian  coordinates  in  a  reference  frame  fixed  in  the  panel,  t  is  the 
time,  Y  the  phase  velocity  with  respect  to  the  x-axis,  o  the  angle  of  propa¬ 
gation  (of  the  wave  normal)  with  respect  to  the  x-axis,  and  k  the  wave 
number-  -  viz.  , 


k  -  Y~“»  (2.2) 

*1 

\j  being  the  wavelength  of  the  surface  wave.  The  amplitude  of  the  wave, 
although  for  convenience  posed  as  unity  in  Eq.  (2.1),  is  assumed  to  be  small 
compared  with  all  other  characteristic  lengths-  -  in  particular,  both  ,he  wave¬ 
length  and  the  boundary  layer  thickness.  We  remark  that  the  assumption  with 
respect  to  boundary  layer  thickness  might  not  be  satisfied  in  many  cases  of 
practical  interest;  but,  of  course,  the  effect  of  the  boundary  layer  on  panel 
flutter  almost  certainly  would  be  unimportant  if  the  amplitude  were  of  the 
same  order  as  or  large  compared  with  the  boundary  layer  thickness. 

Our  model  of  the  boundary  layer  is  an  inviscid,  non-heat-conducting, 
parallel,  shear  flow  that,  in  the  absence  of  the  traveling  wave  disturbance, 
is  characterized  by  the  velocity  profile  U(z),  the  temperature  profile  T(z), 
and  the  free  stream  pressure  pj  (which  is  constant  through  the  boundary 
layer).  The  sonic  velocity  a(z)  and  density  p(z)  then  are  related  to  the 
temperature  according  to 

a2(z)  pj  T  (z) 
a^  p(z)  Tj 


Cf.  Eq.  (3.13),  Ref.  2,  x..  and  y^  therein  being  replaced  here  by  x  and  y. 
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by  virtue  of  the  constant  pressure:  the  subscript  1  on  the  fluid  variables 
denotes  the  free  stream  conditions  outside  the  boundary  layer.  We  neglect 
the  x  and  y  variations  of  all  of  these  quantities  with  the  implicit  assumption 
that  local  values  may  be  used  (as  in  the  boundary  layer  stability  problem). 

If  the  boundary  layer  is  laminar  the  essential  approximation  in  our  shear 
flow  model  is  that  the  viscous  forces  associated  with  the  traveling  wave  are 
negligible  compared  with  the  corresponding  inertia  forces  (the  neglect  of  heat 
conduction  is  almost  certainly  of  less  importance  than  the  neglect  of  viscosity). 
This  approximation  is  similar  to  that  adopted  in  the  inviscid  approach  to  the 
boundary  layer  stability  problem,  although  in  the  prenent  problem  the  approxi¬ 
mation  remains  valid  in  the  neighborhood  of  the  boundary  by  virtue  of  the  finite 
acceleration  there.  As  in  the  boundary  layer  problem,  the  approximation  breaks 
down  in  an  inner  viscous  layer  (Ref.  3,  page  136).  where  the  inertia  forces  tend 
to  sero.  The  result  is  a  singularity  in  the  inviscid  equations  of  motion  (see 
Sections  4  and  5),  which  may  be  said  to  represent  the  inner  viscous  layer  in  the 
limit  of  zero  viscosity  .  We  infer  from  the  known  results  for  the  boundary 
layer  stability  problem  that  this  representation  is  significant  for  small  but 
finite  values  of  the  viscosity  if  the  perturbation  motion  is  unstable. 

If  the  boundary  layer  is  turbulent  and  we  interpret  U.  T,  p ,  a.  and  pj  as 

mean  values,  the  approximations  implied  by  the  model  are  more  severe.  First, 

a 

we  neglect  the  perturbation  Reynolds  stresses  associated  with  first  order 
coupling  between  the  perturbation  flow  and  fluctuations  in  the  original  flow. 
Secondly,  we  implicitly  neglect  cross -correlations  between  the  fluctuations  of 
unlike  quantities  --e.  g., density  and  velocity.  We  suggest  that  the  latter 
approximation  is  not  likely  to  introduce  uncertainties  greater  than  those 
inherent  in  presently  available  profile  data  for  compressible  boundary  layers: 
this  conjecture  is  supported  by  the  observation  that  these  cross -correlations 
would  not  appear  in  the  corresponding  incompressible  problem  and  probably 
represent  only  minor  compressibility  effects  compared  with  the  variation  of 
mean  Alach  number  (which  is  incorporated  in  our  model). 

We  next  remark  that  perturbations  having  the  phase  velocity  V  with 
respect  to  the  steady  shear  flow  will  appear  as  steady  disturbances  to  an 
observer  moving  with  this  same  velocity  V.  with  respect  to  whom  the 
original  flow  will  have  the  velocity  distribution  U(s)-V.  This  allows  us  to 

Gf.  Appendix,  Ref.  7. 
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A 

use  Ward's  result  that  the  perturbation  pressure  p-p^  and  the  perturbation 
velocity  components  v  and  w  (along  y  and  z,  respectively)  may  be  derived 
from  a  potential  $  according  to 


Pi-P  =  «v 

(2.4a) 

P  ( U  -  V)  v  =  *y, 

(2.4b) 

P(U-  V)  w= 

(2.4c) 

where  4  satisfies 

(1  -  M2)  6  +4.  +  <p  -  f|»)(M  4>  +  M  «j>  )  -  0.  (2.5) 

'  '  rxx  yy  Tzz  l  My'  y  Ty  zTz 

with  M  as  the  apparent  (to  an  observer  moving  with  the  wave  at  elevation  z) 
Mach  number 


M  =  M(z)  =  .  (2.6) 

We  emphasize  that  4  is  not  a  velocity  potential  in  the  classical  sense  and  that 
the  flow  is  not  irrotational. 

The  linearized  boundary  condition  to  be  imposed  on  4  at  the  wall  (z  =  zq) 

4f 

is,  from  Eqs.  (2.1)  and  (2.4c), 


We  note  that  if  M(z)  is  introduced  in  place  of  z  as  an  independent  variable 
the  boundary  condition  (2.7)  may  be  regarded  as  imposed  along  the  streamline 
on  which  M=  M.  =  M(z  )  in  the  absence  of  the  disturbance,  after  which  the 
variable  z  may  oe  interpreted  as  a  parametric  streamline  coordinate,  rather 
than  the  linear  distance  from  the  wall.  ^The  validity  of  this  interpretation  may 
be  established  in  two  dimensions,  say  x  and  z,  by  introducing  the  independent 
variables  x  and  4*,  where  4*  is  s  stream  function  for  the  compressible  flow,  in 
place  of  x  and  z  (von  Mises  transformation)  prior  to  the  linearization  of  the 
equations  of  motion.j  This  stratagem  avoids  the  restriction  (that  would  be 
implied  if  z  were  interpreted  directiy  as  a  linear  measure)  to  amplitudes 
that  are  small  compared  with  any  distance  in  the  boundary  layer  over  which 
appreciable  changes  in  velocity  occur;  such  a  restriction  would  be  highly 
undesirable  because  of  the  large  velocity  gradient  near  the  wall.  We  also  note 
that,  while  U  =  0  at  the  wall,  it  might  be  desirable  in  some  applications  to 
impose  the  boundary  condition  (2.7)  at  some  other  streamline- -e.  g. ,  the  edge 
of  the  laminar  sublayer;  this  would  be  permissible  insofar  as  the  distance 
between  this  streamline  and  the  wall  were  negligible  compared  with  the  wave¬ 
length.  hi  effect,  we  assume  that  the  sublayer  moves  with  the  boundary. 
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*  =  p(U-V)^i  =-  ikcoi  a  P(U-V)2i,  z«*.  (2.7) 

Z  JJt  o 

This  boundary  condition  (together  with  the  anticipated  form  of  the  expression 
for  the  pressure  at  s  *  zQ)  suggests  that  we  assume  a  solution  to  Eq.  (2.5) 
in  the  form 

$(x-  Vt,  y,  z)  =  i  Pj  a  2  sec  a  f(z)  t  (x-  Vt,  y)  .  (2.8) 

Substituting  this  in  Eqs.  (2.5)  and  (2.7)  and  introducing 

m  =M(z)  cos  a  =  ^zfi'ecV  ’  (2*9) 

we  obtain 

-  k2(l  -m2)f  =  0  (2.10) 

and 

=  -km2(zQ)  .  (2.11) 

z=zQ 

The  differential  equation  satisfied  by  f(z),  to  wh.ch  the  perturbation 

5 

presE  -re  is  proportional,  was  obtained  by  JLighthill  in  his  studies  on  shock 
wave,  boundary  layer  interaction;  it  is  appreciably  simpler  than  the  differ¬ 
ential  equation  satisfied  by  w  (the  latter  usually  has  been  used  in  boundary 

layer  stability  studies^).  We  also  note  that  if  m  =  0  (U,  =  V)  at,  say,  z  *  z 

•  c 

in  the  boundary  layer,  the  differential  equation  (2.10)  has  a  regular  singularity 
* 

there.  The  immediate  (z)  neighborhood  of  this  singularity  constitutes  the 
inner  viscous  layer,  where  the  inertial  forces,  being  proportional  to  (U  -  V)2, 
cannot  dominate  the  viscous  forces;  however,  the  thickness  of  this  layer 
approaches  aero  for  large  Reynolds  numbers  (Ref.  3,  page  136)  if  the  disturb¬ 
ance  is  self-excited,  and  we  therefore  neglect  it  in  the  following  development 
(except  ineofar  ae  it  ie  represented  by  the  singularity  at  z  =  zc), 

T" -  ■■  "  " 

This  singularity  is  only  apparent  if  (U'/T)'  vanishes  at  z=z  ;  cf.  Eq.  (4.7). 


df 
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The  differential  equation  satisified  by  f(z)  outside  of  the  boundary  layer, 
where  m  assumes  the  constant  value  m^,  reduces  to 


d*2 


(2.  12) 


where 


P 


We  then  have 


...  2.1/2  . 

k(l  -  m  )  =  k 


1  - 


°i-v 
sec  o 


1/2 


(2.13) 


-p  z 

f(z)w  e  ,  *  >  , 


(2.  14) 


where  z  =  Zj  denotes  the  edge  of  the  boundary  layer;  from  this  we  deduce 
the  boundary  condition 


+  P*  =  0,  z  »  Zj  .  (2.  15) 

The  finiteness  condition  that  f  be  bounded  as  z— ♦  oo  requires  Rf  P-  0, 
while  the  radiation  condition  requires  Im  P  to  have  the  same  sign  as  V  -  U  if 
R  IP  =  0.  These  requirements  will  be  satisfied  everywhere  in  the  complex 
V-plane  if  we  draw  the  cuts  for  P  from  the  branch  points  V  =  Uj  ±a^  sec  a  to 
±  cd  and  exclude  points  on  their  top  sides.  Of  course,  is  not  necessary  to 
satisfy  the  finiteness  and  radiation  conditions  everywhere  in  the  V-plane  but 
only  at  those  points  for  which  solutions  to  Eq.  (2.  12)  are  to  be  considered  (thus, 
in  reference  1,  the  cuts  were  drawn  into  Im  V?0  and  only  points  in  Im  V^O 
considered);  in  the  final  analysis,  we  may  deform  the  cuts  in  any  way  that 
ensures  the  existence  of  the  solution  to  the  differential  equation  and  the  satis¬ 
faction  of  the  boundary  conditions.  We  remark,  in  this  connection,  that  a 
periodic  disturbance  of  the  type  (2.  1)  cannot  exist  by  itself,  and  in  a  complete 
initial  value  problem  we  would  be  led  to  consider  an  integral  superposition  of 
such  disturbances  over  a  prescribed  path  in  the  V-plane. 


*In  Section  5,  (3  is  redefined  as  K  (1  *  m j  )  where  K-  =  k6,  and  6  denotes 

boundary  layer  thickness. 
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A  *  sec  a  .  (2.18) 

Our  boundary  value  problem  now  may  be  posed  as:  find  a  solution  to 
the  differential  equation  (2.10),  subject  to  the  boundary  conditions  (2.11)  and 
(2.15),  and  evaluate  it  at  a  =  aQ  . 
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3.  EQUATIONS  OF  MOTION  -  PANEL 

We  present  in  this  section  a  brief  derivation  of  the  eigenvalue  equation 
for  the  panel.  A  more  detailed  derivation  hae  been  given  in  reference  2. 

We  assume  the  existence  of  an  operator  L,  such  that  L(  is  the 
reaction  force  per  unit  area  associated  with  a  transverse  displacement  £  . 

L  may  imply  both  differential  and  integral  operations  with  respect  to  the 
space  coordinates,  but  not  with  respect  to  the  time-  -i.  e. ,  L  £  is  in 
phase  with  £  .  The  equation  governing  the  panel  motion  then  reads 

+  T  *  Pi  ‘  V  (3.1) 


where  p^  denotes  the  perturbation  pressure  of  the  fluid  beneath  (inside) 
the  panel,  p  the  perturbation  pressure  outside  the  panel,  and  <r  the 

A 

panel  mass  per  unit  area.  If  the  compressibility  of  the  internal  fluid  is 
neglected,  p.  represents  a  virtual  inertia  force,  which,  for  the  traveling 
wave  defined  by  Eqs.  (2.1)  and  (2.2),  is  given  by 

pt  =  pjk  (V  cosa)2£,  (3.2) 


where  p.  is  the  fluid  density.  (This  follows  directly  from  the  fact  that 
the  effective  depth  of  an  incompressible  fluid  for  a  surface  wave  motion 
is  1/k;  see  reference  2  for  an  alternative  derivation.) 


The  eigenvalue  equation  corresponding  to  the  traveling  wave  of 

Eq.  (2.1)  may  be  derived  (see  reference  2  for  details)  by  relating  the 

operator  L  to  the  wave  speed  V  that  would  result  in  the  absence  of 

o 

the  aerodynamic  force  (but  including  the  virtual  inertia  force  of  the  internal 

fluid);  thus,  substituting  r  and  p.  from  Eqs.  (2.1)  and  (3.2)  in  Eq.  (3.1) 

with  V  =  V  therein,  we  obtain 
o 


Lf 


f  PjA  2 

<r  +  — —  (k  Vq  cos  o)4r  , 

V  kJ 


(3.3) 


We  emphasise  that  Vq  is  an  eigenvalue  of  Eq.  (3.3)  and  is  a  function  of 
k  and  a . 
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Substituting  from  Eq.  (2,16),  and  L(  from  £q.  (3,3)  in  tha 
com  pi  at*  oquatioB  (3.1),  with  V  i  YQ  in  6  and  pj,  aad  dividing  through 
by  tha  coaffidant  of  V* ,  wa  obtain  tha  aiganvalua  aquation  (for  V) 


(3.4) 


whar* 


*1 

♦  rk 


(3.5) 


ia  tha  aerodynamic -maaa  paramatar,  and  F  and  A  are  definad  by 
Cqa.  (2.17)  aad  (2.18). 
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4.  SOME  GENERAL  CONSIDERATIONS  REGARDING  STABILITY 

The  explicit  determination  of  the  stability  boundaries  for  a  given  panel, 
as  represented  in  the  specification  jor  derivation  from  Eq.  (3.3)j  of  VQ 
as  a  function  of  k  and  a  and  the  mass  parameter  p,  and  a  given  boundary 
layer,  as  represented  by  the  specification  of  velocity  and  temperature 
profiles,  requires  the  simultaneous  solution  of  the  differential  equation  (2.10), 
subject  to  the  boundary  conditions  (2.11)  and  (2.15),  and  the  transcendental 
equation  (3.4).*  This  represents  a  formidable  problem,  which  is  especially 
complicated  if  Y  =  U  within  the  boundary  layer,  giving  rise  to  a  singularity 
of  the  differential  equation  (2.10).  In  order  to  clarify  the  role  of  this  singu¬ 
larity  with  respect  to  panel  flutter,  we  develop  in  this  section  an  implicit 
solution  to  the  boundary  value  problem. 

We  first  rewrite  Eq.  (3.4)  in  the  form 


where  the  phase  of  the  radical  is  defined  such  that  V  =  Vq  at  n  =  0.  It 
then  follows,  neglecting  structural  damping  (Vq  real),  that  the  imaginary 
part  of  V  will  be  opposite  in  sign  to  the  imaginary  part  of  F,  so  that 
instability  requires  Im  F>0. 

We  obtain  an  implicit  solution  for  F  by  multiplying  Eq.  (2.10)  through 
by  m~27  (where  the  overbar  denotes  the  complex  conjugate),  integrating 
from  x  =  *0  to  x  =  2^ ,  integrating  the  term  T  (m“  f*)*  by  parts,  imposing 
the  boundary  conditions  (2.11)  and  (2.15)  to  eliminate  f*  at  the  limits,  and 
solving  for  T  (xq);  the  end  result  is 


»  s  r<*o>  - 


f(V 


2  -2  :/1  2 
ml  ^  -  ml 


1/2 


+  i  rat 


f’U)  |2  +  k2  (1  -  m2) 


ds 

2 

m 


(4.2) 


Alternatively,  Eq.  (3.4)  may  be  combined  with  Eq.  (2.11)  to  obtain  a  homo¬ 
geneous  boundary  condition  on  f  at  i  :  ic,  This  formation  leads  to  a  more 
symmetric  boundary  value  problem  but  is  less  satisfactory  for  the  approxi¬ 
mate  solutions  treated  herein. 
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We  will  aaeume  that  tha  phaae  angle  of  V  ia  amall  ( V  approximately  real): 
then  m  will  he  approximately  real,  and  the  imaginary  part  a  of  the  firat  and 
aecond  tarma  on  the  right-hand  aide  of  Eq.  (4.2)  muat  be  derived  chiefly  from 
the  radical  (1  -  m|)^  and  the  aingularity  at  m  *  0,  reapectively. 

2  1/2 

The  phaae  of  (I  -  )  '  *  0/k  for  V  approximately  real  ia,  by 

definition,  0  for  |mj|  <  1  or  4  t/2  for  *  >  1 ,  correaponding  to  our 

previoua  choice  of  branch  cuta  for  0  jaee  remarka  following  Eq.  (2.1 5]j  . 
Then,  writing 


>  1 


(4.3) 


we  infer  that  the  contribution  to  Im  F  ( s  -ImT)  of  the  firat  term  in 
Eq.  (4.2)  ia  F%,  0.  or  -  Fft  for  V<  Uj  -  A,  Ut  -  A<V<Uj  +  A.  or 
V>(J|  +  A,  reapectively,  if  V  liea  doae  to  the  real  axia  but  not  too 
clone  to  Uj  4  A.  The  integral  in  Eq.  (4.2)  will  contribute  an  imaginary 
part,  aay  Fb,  to  F  only  if  0<V<Uj,  in  which  caae  m  will  vaniah  at 
the  point  »  ■  «c>  It  followe  that 

ImF*F%+Fb,  0<V<Ui-A 


*  Fb,  UrA<V<Ut 

*  0,  U|  <  V<  Uj  +  A 


(4.4) 


i  -Fa.  Uj  +  A<V. 


We  remark  that  if  V  done  approximate  Uj  *  A,  aa  in  Section  6,  the 
contribution  of  the  firat  term  in  Eq.  (4.2)  to  ImF  muat  be  modified, 
but  thia  doea  not  affect  the  contribution  of  the  aecond  term,  which  remaina 
aa  Fb  or  0. 

The  nature  of  the  negative  damping  repreaented  by  Fft  ia  eaaentially 
the  aa me  aa  that  predicted  in  the  abaence  of  the  boundary  layer  and  diacuaaed 
in  more  detail  in  reference  2.  The  preaence  of  the  boundary  layer  reducea 
thia  component  by  partially  iaolating  the  aurface  wave  on  the  panel  from  the 
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direct  action  of  the  high-speed  airflow  outside  the  boundary  layer  (in  particular, 
the  panel  is  exposed  directly  only  to  a  subsonic  flow,  even  though  the  free- stream 
flow  is  supersonic). 


We  may  deduce  a  simplified  expression  for  by  noting  from  the 

differential  equation  (2.10)  that  f'  (s)  -  O(m)  and  f(s)  *  f(s  )  +  0(m2)  as 
*  c 
z-*z  i  accordingly, 

% 


(4.5a) 


=  -k  6 


f(sc)  2  1m  Kt, 


(4.5b) 


where  Kj,  as  defined  by 


1)  da, 


(4.6) 


is  the  integral  appearing  in  the  solution  for  f  j  (a)  discussed  in  Appendix  A. 

The  sign  of  Im  K,  evidently  depends  on  how  the  singularity  is  circumvented; 

this  is  an  extremely  delicate  question,  but  it  has  been  examined  thoroughly 

by  Lin  (Ref.  3,  ch.  8),  and  we  merely  state  his  result  that  the  path  of  inte- 

*% 

gration  must  be  indented  over  the  point  a  =  «c,  whence 


vka 


2 

sec 


d_ 

da 


(4.7) 


where  the  subscript  c  implies  evaluation  at  a  =  ac  (see  Appendix  B  for  details). 

The  integral  Kj  occurred  in  the  laminar  boundary  layer  stability  problem, 
and  Lees  and  Lin^  have  shown  that  if  the  original  flow  absorbs  energy  from  a 
small  disturbance  in  the  boundary  layer  profile 


*  Cf.  Eq.  (5.4b). 

* 

We  note  that  the  sign  of  the  wave  number  usually  assumed  in  boundary  layer 
studies  is  opposite  that  chosen  here,  which  has  the  effect  of  interchanging  i 
and  -i;  thus,  the  path  of  integration  elected  by  Lin  actually  is  indented  under 
the  singular  point  z  =  a  . 


<  0. 
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> 


_d 

da 


■■a. 


(4.8) 


It  doaa  not  follow  that  Eq.  (4.8)  la  a  necessary  and  aufficiaat  condition  for 
tha  atability  of  tha  boundary  layer,  but  it  aaama  likely  jon  the  baala  of 
empirical  conaiderationa,  aa  well  aa  Eq.  (4.8))  that  it  will  be  eatiefied 

J 

for  both  laminar  and  turbulent  boundary  layer  a  in  thoae  aerodynamic 
regimea  for  which  panel  flutter  preaenta  a  serloua  problem.  We  may 
anticipate,  therefore,  that  F^  will  be  poeitive,  and  we  infer  from  thia 
conaideration  that  the  boundary  layer  could  be  deatabiliaing  with  reapect 
to  aurface  wave  a  having  wave  apeeda  in  the  interval  0*  V<  Uj  and 
amptitudea  email  compared  with  both  boundary  layer  thickneaa  and  wave¬ 
length;  however,  the  net  effect  in  the  interval  0<  V<-U.  -  A  ia  more  likely 

*  * 

to  be  atabilising  in  conaequence  of  the  reduction  of  .  The  deatabiliaing 
mechanism  repreeented  by  F^  appeara  to  be  cloaely  related  to  that  atudied 
in  the  boundary  layer  atability  problem  (aee  reference  3,  Section  4.4;  alao 
reference  7). 

A  rough  approximation  to  the  value  of  a  for  a  turbulent  boundary 
layer  may  be  deduced  from  the  well-known  i/7  power  law- -vis. , 


(4.9) 


where  8  ia  the  boundary  layer  thickneaa.  We  anticipate  that  panel  flutter 
will  be  moat  aerioua  when  V  approximate  a  or  ia  leas  than  Uj  -  a}  (aa  in 
the  absence  of  a  boundary  layer;  cf.  reference  2  and  Sections  6  and  7  below). 
Setting  U  *  Uj  -  Sj  and  a  *  *c  in  Eq.  (4.9)  then  yields 


(4.10) 


We  note  that  F  *  0  in  eubsonic  flow,  and  instability  then  would  have  to  be 
charged  entirely  to  the  ahear  flow  in  the  boundary  layer.  It  aeema  likely 
that  IV  would  be  too  small  to  render  practical  aircraft  structures  unstable, 
but  it  may  afford  an  explanation  of  the  generation  of  gravity  wavea  on  a 
liquid  aurface.  ? 
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which  hat  the  value*  5  •  10'*  and  8  *  10*'  for  M,  ■  1.5  and  2,  respectively; 

these  values  are  so  small  that  they  might  place  s  in  the  laminar  sublayer, 

c 

where  the  1/7  power  law  no  longer  would  be  valid,  but  they  may  be  assumed 
to  give  the  correct  order  of  magnitude.  Of  course,  (*c  -  sq)  /  6,  as  given  by 
Eq.  (4,10),  does  not  remain  small  for  large  Mj,  but  it  seems  likely  that 
F.  then  would  be  unimportant  compared  with  F  • 

We  add  that  (sc  -  *Q)  /  6  probably  would  approximate  (1  -  1/Mj)  for  a 
laminar  boundary  layer,  but  a  laminar  layer  almost  certainly  would  be  so 
thin  as  to  have  no  effect  on  panel  flutter. 
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5.  THE  BOUNDARY  VALUE  PROBLEM 

We  now  return  to  the  boundary  value  problem  poeed  by  the  differential 
equation  (2.10)  and  the  aeeociated  boundary  condition*  (2.11)  and  (2.15). 
Introducing  the  dimensionleee  variable 


and  the  boundary  layer  thicknes*  parameter 


(5.1a) 


6  =  -  *0  .  (5.1c) 

in  Eqa.  (2.10),  (2.11),  and  (2.15),  we  obtain 

m2  jm'2  f'j  -  *2(1  -  m2)  f  *  0, 

f0  “  *  *mQ  * 

and 

fg#  +  Pft  s  0,  (5.3b) 

where  the  prime*  imply  differentiation  with  respect  to  t, ,  the  subscript* 

0  and  1  imply  evaluation  at  £  =  0  and  1,  and  (3  is  defined  by  Eq.  (2.13) 
if  k  is  replaced  by  K  therein. 

The  differential  equation  (5.2)  has  a  regular  singularity  at  £  *  £  , 

c 

where  m(£c)  *  0.  The  exponent*  of  this  singularity  are  3  and  0,  and  we 
find  that  the  corresponding  solution*  in  it*  neighborhood  are  of  the  form 

w,iu  .  t  i  -  ;e>3  -  &  u  -  ;c>4  +  ou  -  ic>5 


Of  ;il,  (5.2) 

(5.3a) 


(5.4a) 


( 
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and 


W2U) 


2 

Ji*~  w 

3 


2 

,  (c)iog  [<;-  cc»«‘1,r]+ 1  - ■*-  u- c.c)2  +  <xt-  4ci 


(5. 'lb) 


where 


-m 


q  = 


(5.5) 


m 


The  solution  w^  is  everywhere  real  for  real  values  of  £  -  £c,  while  w2 
has  been  normalised  to  have  the  imaginary  part  0  or  (-tri q/C2/3) Wj  as 


£  <£  or  £  ;>  £c»  respectively  (since  the  logarithmic  branch  point  must  be 
circumvented  by  indenting  over  £  =  £c).  The  Wronskian  of  this  pair  of 
solutions  is 


WiW2 


‘  W2Wi 


3mfc(rj 
- 72  ‘ 


m 


(5.6) 


We  also  find  it  convenient  to  introduce  a  second  pair  of  solutions,  f^  and 
that  satisfy  the  boundary  conditions 


fI0  =  i.  (5.7a) 

f’10  *  0,  (5.7b) 


(5.7c) 


and 


<5.7d) 


where  the  first  and  second  subscripts  identify  the  solution  and  the  point  of 
evaluation  respectively.  The  corresponding  Wronskian  is 


w|fr  f.j  =  m2(£). 


(5.8) 
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Invoking  the  boundary  condition*  (5.7)  on  linear  combination*  of  and  w2 
and  making  use  of  the  Wronsklan  (5.6),  we  obtain 


m 

ft(c)  -  ~r  [So  w2<;) '  w2o  wi  <c>] 


(5.9a) 


'  4 

f2(C)  =  — ---  [w20  Wj  (C)  -  wl0  w2  (t)J  •  (5.9b) 


The  required  solution,  subject  to  the  boundary  condition*  (5.3),  now 
may  be  expressed  in  terms  of  either  and  or  f j  and  We  find 


-  Km, 


2  Kj  +  ft»u>  W2U)  -  (w2t  +  Pw21)  wt  (£) 

°  Ki  +  Pwli)  w20  '  (w2i  +  ^2i>  wi0 


(5.10a) 


f  *  Ffjte)  -Kf2(C), 


where 


F  =  K 


fkl  +  ?f21 

<il  +  *£ll 


(5.10b) 


(5.11) 


We  require  only,  the  imaginary  part  of  F  throughout  much  of  the  sub¬ 
sequent  analysis.  This  can  be  calculated  directly  from  Eq.  (5.11),  but  it 
may  be  advantageous  to  consider  separately  the  components  F  and  F,  , 
as  defined  by  Eqs.  (4.3)  and  (4.5b);  these  depend  only  on  the  absolute  value 
of  f  and  do  not  require  the  calculation  of  phase  (which  might  be  rather 
small  in  some  applications).  Setting  £  *  1  in  Eq.  (5.10b),  simplifying  the 
result  with  the  aid  of  the  Wronskian  (5.8)  evaluated  at  £  s  l,  and  substi¬ 
tuting  in  Eq.  (4.3),  we  obtain 


Fa  * 


K2mf(mf  -1)1/2 


fll  +  pfli 


IHj  >  1  . 


(5.12) 
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Setting  j;  =  £c  in  Eq.  (5.10b),  expressing  the  denominator  of  the  result  in 
terms  of  f^  with  the  aid  of  Eq.  (5.9a),  and  substituting  in  Eq.  (4.5b),  we 
obtain 


S'  -  V3'4 

Fb  "  mc 


w'tl  +  Pw,  1  2 


11 


’I'n^iii 


ImKj . 


(5.13) 


The  solutions  fj,  f^,  and  are  developed  as  power  series  in  K  in 
Appendix  A.  The  first  approximations  to  the  quantities  required  for  the 


evaluation  of  F  and  F.  are 
a  d 


fu  =  1+0(KZ),  t'n  =  #C2mfKt  jl+0(K*)j. 


wn  =  ~rr  ! 1  +  ° ( k2)  .  w‘u  =  ^  [i+o(K2)j, 

m  m_  L  J 


(5.14a,b) 


(5.15a,b) 


where  is  the  integral  introduced  in  Eq.  (4.6)  and  evaluated  in  Appendix  B, 
and  Jj  is  defined  by  Eq.  (Alii)  in  Appendix  A.  Substituting  these  approxima¬ 
tions  in  Eqs.  (5.1?.)  and  (5.13)  yields 


2.2  ,,1/2 

m,  (m,  -  1) 


-T-m-1  -r-rz  [*  +  °<‘2>j 

i(m2  -  1)1/2  +  Kmf  K{  \£  1  J 


=  0  , 


and 


F. 


I  2  x  I/M  2,  i/2  _ 

'irij  +  K(1  -  m1 )  Jj 

j  (1  -  m[) 172  +  Km2 


-Kmj  ImKj 


Im  K, 


1  +  0(K  ) 


mt>l,  (5.16a) 


m^  1 ,  (5.1 6b) 


(5.17a) 


(1  -  m2)  1/2  +  K  m2 Kj 


i  +  o(K2.  K,  1  -  m2  ) 


(5.17b) 
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The  approximation  (5. 17b)  anticipates  (on  the  basis  of  the  reeulte  for  =  0; 

•ee  reference  2)  that  supersonic  panel  flutter  should  prove  most  serious  for 
values  of  m j  only  slightly  larger  than  unity  for  sufficiently  small  values 
of  K. 

The  approximations  to  and  provided  by  Eqs.  (5. 16)  and  (5. 17b), 
together  with  the  values  of  Kj  given  by  Figs,  la  and  lb  (for  a  =  0  and  an 
insulated  boundary)  are  plotted  in  Figs.  2a  through  2c  and  Figures  3a  through  3d, 
for  a  0.1,  1.2,  1.5,  and  2.  0  and  K  s  0.  5,  1.0,  and  2.  0.  Also  plotted 
(labelled  "exact")  are  some  results  that  were  obtained  by  numerical  integration 
of  the  differential  equation  (5.  2)  on  a  digital  computer.*  The  agreement  doubtless 
could  be  improved  by  including  terms  of  higher  order  in  K  in  the  approximate 
results  for  F&  and  F^  or,  in  some  cases,  merely  by  using  Eq.  (5.  17a)  in 
place  of  Eq.  (5.  17b),  but  the  errors  in  the  approximations  (5. 16a)  and  (5. 17b) 
are  not  likely  to  exceed  those  already  present  in  our  simplified  model  of  the 
boundary  layer. 

We  emphasise  that  the  results  (both  "exact"  and  "approximate")  for 
are  not  valid  for  V/Uj  less  than  about  0.  3  in  consequence  of  the  breakdown 
of  the  one-seventh  power  law  assumed  in  the  calculations.  An  adequate 
correction  factor  for  F^  could  be  achieved  by  forming  the  ratio  of  Im 
for  the  actual  U  (s)  and  T  (s)  to  that  based  on  the  one- seventh  power  law, 
using  Eq.  (B8c).  The  results  for  Fft  do  not  appear  to  be  very  sensitive  to 
changes  in  the  profile. 


The  procedure  was  to  obtain  first  the  regular  solution  w j  by  integrating  away 
from  the  singularity,  starting  the  solution  by  power  series  expansions  near 
(  *  L  and  then  continuing  by  the  Rungl-Kutta  method  over  the  range  (2“*  „  1). 
The  singular  solution  then  was  expressed  as 

w2  =  (qK2^)  wjlog  jjfc -&c)e”l1j  +g(£), 

the  inhomogeneous  differential  equation  for  g  determined  by  substitution  in 
the  original  differential  equation,  and  the  regular  function  g  obtained  by  the 
procedure  outlined  for  wj.  In  those  cases  where  (  was  outside  of  the  range 
of  integration,  f  j  was  determined  directly  by  integrating  out  from  (  =  2“  *  . 
The  details  and  programming  were  worked  out  by  Dr.  Samuel  Conte  and 
Mr.  David  Bussard;  the  computer  was  a  Remington-Rand,  Univac  Scientific 
Model  1103A. 


as  (B16)  -  (B20) 
0  ths  Niultl 
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v/U, 


Figure  2b.  F^  versus  V/Uj  for  Mj  =  1.5. 


GM-TR-299 

Pag*  26 


.5. 


6.  APPROXIMATE  STABILITY  ANALYSIS 


We  consider  in  this  section  the  approximate  calculation  of  the  negative 
damping  ratio  (in  the  notation  of  reference  2) 

=  -  Lm  V/Ri  V  (6.1) 

on  the  assumption  that  the  second  term  in  the  braces  of  Eq.  (4.1)  is  sufficiently 
small  compared  with  unity  to  permit  the  approximation 

f 

V=Vo{‘- 

Eq.  (6.2)  may  be  regarded  as  the  first  approximation  based  on  a  power  series 
expansion  in  the  parameter  p,  and  it  will  be  valid  for  p  <<  1  insofar  as  F 
remains  bounded.  Now  p  [  see  Eq.  (3.5)]  almost  certainly  would  be  small 
in  practical  applications,  while  |  F  |  can  be  large  only  if  both  J  V  -  (Uj  ±  A)| 
and  K  «  1;  the  latter  contingency  is  examined  in  the  following  section,  but  we 
remark  here  that  it  is  apt  to  be  important  only  when  K  is  of  the  same  order  of 
magnitude  as  p1^. 

Substituting  Eq.  (6,2)  in  Eq.  (6.1)  and  evaluating  Im  F  from  Eq.  (4.4) 
yields 

6o*  l*  N7J <F*+Fb>-  °<VU1  -  A  '6-3> 

in  the  region  of  principal  interest.  The  maximum  value  of  6q  with  respect 

to  V  almost  certainly  will  occur  in  this  region  if  (V  )  <  U,  -  A,  but 

o  7  *  '  o'min  1  ’ 

F  =0ifU.-A<V  <  U, ;  if  V  >U,  panel  flutter  of  the  type  considered 
a  x  o  x  o  i 

here  is  not  possible. 

Perhaps  the  most  expedient  approach  to  the  determination  of  the  maximum 
value  of  with  respect  to  variations  of  K  and  a  is  to  adopt  the  approxima¬ 
tions  of  Eqs,  (5,16)  and  (5.17b)  for  Fa  and  F^,  determine  the  corresponding 

value  of  (6  )  ,  and  then  estimate  the  approximate  corrections  to  F  and  F, 

o  max  "  a  b 

from  the  results  of  Fig.  2.  Substituting  Eqs.  (5.16)  and  (5,17b)  in  Eq.  (6.3) 
yields 
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1  [  A  \ 

=  2  a  hr~  I 

V  o  j 


m  2  Im 


1 


1  +  Km  2  Kj 


(6.4a) 


1  /U1  ^2 

i 11  nr  " 1 


m  2  -  1  -  Km  2  Im  Kj 


. -  ' i - J — LJ _ 

^n2  -  1  -  Km  2  Im  Kx  j  +  Km  2  R  i  Kj 


j-g — •  1 ,  (6.4b) 


We  recall  that  [Eq.  (2.9)] 


U1‘ Vo 

ml  =  I  Tf2  )  C0S  a  =  M1 


V  A 

1-eA  c 


7JT J  COB  a  * 


(6.5) 


while  Kj  depends  on  V^/Uj,  Mj,  and  cos  a;  as  suggested  in  Appendix  B, 
the  error  in  assuming  that  Kj  depends  only  on  Mj  cos  a,  rather  than  M ^ 
and  cos  a  separately,  usually  will  be  small  (the  essential  approximation  is 
the  replacement  of  Mj  by  Mj  cos  a  in  the  temperature  profile). 

We  consider  as  a  more  specific  example  a  pressurized  cylindrical  shell 

for  which  the  internal  pressure  is  sufficiently  high  to  render  negligible  the 

effects  of  bending  on  VQ,  which  then  is  given  by  [  Ref.  2,  Eq.  (5.5)  with 

D  =  0  and  N  =  2N  a  p.R  I 
y  x  J 


V  2  =  — 
o  cr 


Eh 


cos  a 
kR 


+  p^R  j  j  +  tan2  a 


We  introduce  as  the  independent  variables 

x  =  VQ/Uj,  y  =  cos' 


(6.6) 


(6.7a,  b) 


and,  as  a  measure  of  the  wavelength, 

z  =  Zj  /kR,  Zj  =  ^2  Eh /p.R 

Substituting  (6.7a,  b)  and  (6,8a,  b)  in  (6.6),  we  obtain 


(6.8a,  b) 


z=  ,i 
y 


i  +! 


m 


2 

“2 

y 


(6.9) 


where 
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m 


(6.10) 


denotes  the  minimum  value  of  x  (x  =  x  at  a  =  0  and  k  =  oo).  The  boundary 

m  r 

layer  thickness  and  aerodynamic  mass  parameters  I  the  latter  given  by 

-i 

Eq.  (3.5)  with  f.  *0  therein  jthen  may  be  expressed  as 

VC  *  k6  =  tl  (6/R)  a'1  (6.11) 

and 

p  a  pl  /<rk  =  (pjR/traj)*  .  (6.12) 

Noting  also  that 

m  2«  M  2  (1  -x)2  y  (6.13) 

and,  on  the  basis  of  the  aforementioned  assumption  for  Kj  , 

m2  Kj  *  m  2  K1(Vo/U1.  Mj  cos  a)=m2  Kj  (x,  Mjy1/2)  ,  (6.14) 


Eq.  (6.4b)  may  be  transformed  to 

_  ,  1 

“  1  jMl  ^  "  X)2V2  "  1  +*j(5/R)*"1  m2Kj(x,  Mjy1  ^2) 

(6.15) 

We  require  the  maximum  value  of  5q,  as  given  by  Eq.  (6.15)  in  con¬ 
junction  with  Sq.  (6.9).  with  respect  to  independent  variations  of  x  and  y 
over  the  intervals 


x_  4  x  £  1 
m 


1  +  (x/x  )2 
m 


and 


6  y  £  1  . 


(6.14a,  b) 
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The  required  calculation*  would  be  difficult  to  carry  out  analytically,  but  they 
are  quite  simple  on  a  high-speed  computer. 

If  the  boundary  layer  is  relatively  thin  (K«  1)  we  find  that  (6Q)max  i* 
likely  to  occur  for  values  of  x  and  y  such  that  exceeds  unity  by  only 

a  small  amount  and  the  two  terms  in  the  braces  of  Fq.  (6.15)  are  approxi¬ 
mately  equal  in  magnitude.  The  imaginary  part  of  Kj--and,  therefore,  the 
negative  damping  effect  of  profile  c.urvature--then  is  likely  to  prove  rather 

unimportant,  although  its  real  part  plays  an  essential  role  in  preventing  the 

* 

infinity  that  would  be  indicated  by  Eq.  (6.15)  at  nrij  =  1  if  K— 0  (in  which 
case  the  analysis  of  the  following  section  would  be  required).  The  approxi¬ 
mations  of  Eqs.  (5,16)  and  (5.17b)  evidently  are  entirely  adequate  for  this 
case. 

If  the  boundary  layer  is  relatively  thick  (say  K  >  1  or  2)  and  xm  is  not 
too  small,  it  appears  that  (^Q)max  i®  likely  to  occur  at  nij  =  1  and  o  =  0. 

The  instability  in  this  case  would  be  associated  entirely  with  profile  curvature, 
and  it  therefore  would  be  important  to  improve  the  approximation  to  F^, 
relative  to  that  of  Eq.  (5,!7b),  and  to  investigate  the  effects  of  possible 
departures  from  the  one-seventh  power  law. 

Calculations  based  on  the  results  of  this  section  have  been  carried  out 
for  pressurized,  monocoque  shells  filled  with  either  gas  or  liquid  and 
compared  with  calculations  based  on  the  results  of  reference  2„  It  appears 
that  the  degree  of  instability  for  supersonic  flutter  may  be  reduced  by  an 
order  of  magnitude  (roughly  a  factor  of  5  -  10)  in  consequence  of  boundary 
layers  for  which  IO  l/2. 


The  effect  of  the  boundary  layer  then  may  be  compared  with  that  of  small 
damping  in  limiting  the  response  of  a  simple  oscillator  at  resonance. 
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7.  STABILITY  ANALYSIS— VERY  THIN  BOUNDARY  LAYER 


The  anal y el •  of  tha  preceding  section  breaks  down  in  the  neighborhood 
of  »j  «  I  if  K  is  of  the  same  order  as  p1^3  (of  course,  if  K  is 
sufficiently  small  compared  with  p  the  effects  of  the  boundary  layer  on 
panel  flutter  will  be  negligible,  but  turbulent  boundary  layers  are  not  apt  to 
be  this  thin  at  supersonic  speeds).  The  assumption  that  both  K  and  jm^  -  l] 
are  small  allows  us  to  approximate  F  by  £•••  Eqs.  (5.1  i)  and  (A8a)»(A9b); 
also  cf.  Eqs.  (5.16)  and  (5.17a),  with  which  this  approximation  is  consistent] 


F 


T7T 


(1 


m2) 


KmjK 


3 


(7.1) 


Substituting  F  from  (7.1)  in  (3.4)  and  eliminating  through  (2.9)  and 
(2.18).  we  obtain,  after  some  algebraic  manipulation. 


where  the  dependence  of  Kj  on  both  U/A  and  V/U^  has  been  explicitly 
denoted. 

We  next  introduce  the  dimensionless  variables  v  and  £  ,  representing 
the  departures  of  V  and  VQ  from  Uj  -  A,  according  to  (cf.  reference  2) 

V  «  (Uj -A)  [l+(l/2)v(p/2)2/3(Mo- 1)"5/3]  .  -3w/2<arg  v  -  w/2  (7.3a) 

and 

Vo  *  (U1  -  [l  +  <3/2)6  <p/2)2/3  <Mc  -  l)'5/3j  •  (7.3b) 


where  the  permissible  range  of  arg  v  is  deduced  from  the  restrictions  on 
V  and  0  |Eq.(2.13)],  £  is  real  (structural  damping  being  neglected),  and 
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Mq  denotes  the  (free  stream)  Mach  number  relative  to  an  observer  moving 
with  the  wave  front- -viz.  , 

M  =  U./A  s  U,  cos  o/a  .  (7.4) 

a  *  1  co 

We  also  find  it  convenient  to  introduce  the  boundary  layer  parameters  6 
and  n  according  to 


and 


e 


1/3 

*1  (Ma  . 


n 


8  . 


(7.5a) 


(7.5b) 


We  will  use  the  parameter  6  initially,  but  we  find  tj  more  convenient 
in  the  end  results.  Substituting  Eqs.  (7.3a,b)  and  (7.5a)  in  Eq.  (7.2)  and 

neglecting  higher  powers  of  p  then  yields 

(v  -  30  <v1/2  +  30)  +  2  =  0.  (7.6) 

We  observe  that  the  importance  of  the  boundary  layer  in  this  approximation 
is  determined  by  the  size  of  3|0  (relative  to  v  *^2  ;  in  particular,  K<<  1 
is  not  a  sufficient  condition  for  the  neglect  of  the  boundary  layer. 

Rewriting  Eq.  (7.6)  in  the  form 

(  vl/Z  +9  )3  -  3  AB(  vi/Z  +  0)  +  (A3  +  B3)  =  0,  (7.7) 


where 


A 

B 


=  j(i  -  3€e  +  e3)  *  [d  -  3  $e  +  e3)2  -  (&  +  e2)j  1/2J 


1/3 


(7.8) 


we  have 


v 


1/2 


» 


+  9  =  A«  +  B?  , 


(7.9) 
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where 4 the  three  roots  of  the  cubic  equation  (7.7)  correspond  to  the  three  cube 
roots  (  «  )  of  >1  (7  is  the  complex  conjugate  of  «).  This  result  is  valid  for 
complex  values  of  both  (  and  9 .  but  the  explicit  determination  of  the 
imaginary  part  of  v ,  on  which  the  damping  ratio  depends  (assuming  p  to 
be  approximately  real),  is  possible  only  if  both  £  and  9  are  assumed  to 
be  real;  on  this  last  assumption  we  find 

Imv  =  -  (£>2)  |A  -  B|  •  |A+  B  -  29  |  (7.10a) 

(1  -  3*9  +  93)  *  (£  +  93)2 


=  0.  (7.10b) 

If  9  is  not  real  1m  v  must  be  determined  numerically  from  Eq.  (7.9). 

The  maximum  value  of  -Im  v  with  respect  to  the  spectral  parameter  £ 

is  determined  in  Appendix  C.  We  designate  the  ratio  of  this  maximum  to 

its  value  in  the  absence  at  the  boundary  layer  (3 l^2/21^3  at  9  =  0  and 

£  *  0)  as  h(q),  which  is  given  by  Eqs.  (C7)  and  (C3)  and  is  plotted  in  Fig.  4; 

t  ,  the  value  of  £  at  which  this  maximum  occurs,  is  plotted  in  Fig.  5.  The 
m  p 

maximum  (with  respect  to  £  )  negative  damping  ratio  then  is  given  by  Icf. 

Eq.  (4.13),  reference  2^ 

6  *  (-Im  V/RiV).  .  *  (31/2/4)p2/3(M  -  if 5/3  h(  ,j ). 

°  5  *m  Q  (7.11) 

We  emphasise  that  the  results  for  r\  0  may  not  be  directly  significant,  since 
Kj  and  r\  will  be  positive  real  if  the  singularity  at  in  is  excluded,  while 
if  it  is  included  Kj  and  q  must  be  complex;  however,  these  results  may 
be  useful  in  determining  a  first  approximation  to  6Q  for  complex  tj. 

2 

We  find,  by  numerical  comparison,  that  the  approximations  £m  *  -  r\ 
and  h  *  (2*^ 3/3tj)  are  adequate  for  tj > l  [Eqs.  (CIO)  and  (Cll)j .  Substituting 
these  approximations  in  Eqs.  (7.3b)  and(7.11)  and  eliminating  n  via  Eq.  (7.5) 
yields 

6q  *  ,j./4  KKj  (Mq  -  l)2  . 


■1  >  1 


(7.12) 
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The  real  and  imaginary  parts  of  K( ,  as  evaluated  in  Appendix  B,  are 

4  * 

plotted  in  Fig.  6  for  V/Uj  si.  (cf.  Fig.  1).  As  suggested  above, 

the  error  in  replacing  Mj  by  cos  o  for  a  ji  0  is  not  likely  to  be  large. 


IN 
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APPENDIX  A 

SOLUTION  OF  DIFFERENTIAL  EQUATION 


We  require  solution!  to  Eqs.  (5.  7)  and  (5.  8)— vie. , 


2  d  | 

m  “3T  1 

/■ - - 

3 

a 

IV 

-  f  .  0, 

(Al) 

fio  =  lr 

fio  =  °- 

(A2a) 

o 

ii 

o 

fU 

Vi 

f20  *  m0' 

(A2b) 

where  the  second  subscript  indicates  the  value  of  £. 

A  solution  that  ie  useful  for  small  or  moderate  values  of  K  may  be 
obtained  by  introducing  the  expansions 


Lii.K2)  =  2°  K2nf.<n)(C),  i  =  1,2 

1  „  -n  1 


n=0 


(A3) 


in  Eq.  (Al),  which  then  yields 


d 

dl 


f 


m 


,  df. 

■2  l 


(n) 


=  (m-2  -  1)  f. 


(A4) 


We  may  satisfy  the  boundary  conditions  of  Eq.  (A2)  with  the  zero'th 
approximations 

f,'0'  =  i 


(A5a) 


and 


n* 

(0)  _  I 

*  “  1 
JO 


rr>2dC. 


(A5b) 
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The  succeeding  terms  (nel)  in  Eq.  (A3)  then  must  satisfy  the  homogeneous 
boundary  conditions 

f.<n)(0)  =  r<n)'(0)  *  0.  (A6) 

Integrating  Eq.  (A4)  subject  to  these  conditions,  we  obtain  the  recursion 
formula 

ft  f* 

f.(n)  =  I  m2d;  '  (m“2- 1)  f.tn“1>dt,  nil.  (A7) 

Jo  jo 

We  require,  in  the  end  [cf.  Eq.  (2.  24)J],  the  values  of  f . ,  f.  ,  f-,  and 

1  t>  8  i  l 

f ,  only  at  (  =  1.  Following  Lees  and  Lin  *  (but  with  slightly  revised 
notations),  we  may  express  these  in  a  series  of  definite  integrals  according  to 


£n 

-  Z  K2nK 

-  K2n- 

(A8a) 

n=0 

fi'i 

-  m2  T  l< 2n+2  Y 

-  tn  S  K  K2  j  , 

(A8b) 

n=0 

=  m  K  H2n, 

n=0 

(A9a) 

(21 

-  K2"H2„+1  . 

(A9b) 

Ko  =  l* 

P1 

Kj  =  i  (m_2-l)  dC  . 

r1  i* 

K2  =  I  m2dt  I  (m‘2-l)di, 

Jo  Jo 
n  ft 

K3  =  j  (m‘2-l)d;  I  m2d;  (m‘2-l)dC . 

dO  do  JO 


(A10) 
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I 


(All) 


The  numerical  evaluation  of  these  integrals  for  a  laminar  boundary  layer  has 
been  discussed  in  some  detail  by  Lees.  Appendix  B  contains  a  discussion  of 
the  evaluation  of  Kj  (the  most  important  of  these  integrals)  for  a  turbulent 
boundary  layer. 


The  analytic  solution  ,  as  introduced  in  Section  5,  also  may  be  expanded 
according  to  Eq.  (A3).  Integrating  Eq.  (A4)  subject  to  the  boundary  conditions 


w. 


?  <t  *  tc)' 


(A12) 


we  obtain 


2 

m  dC , 


(A13a) 


ri  ,  /*l 

a  /  m  d;  I  (m-2 -1)  Wj  ^n“  ^dt. 

K  K 


(A13b) 


We  require  w^(l)  and  w^(l),  which  may  be  placed  in  the  forms 


w 


11 


m'c 


2  ?L/2nT 

ml  =  K  J2n ' 


( A14a) 


w, .  *  — , 

11  m'  ‘ 
c 


00  •> 

E  Krnj 


2n+l  ’ 


(A14b) 
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APPENDIX  B 
EVALUATION  OF  Kj 

We  consider  the  integral 


6 

Kj  =  6 mlJ  (m"2  -  1 )  d*  ,  m  (*) 


U(«)  -  V 

a(z)  Bee  a 


Introducing  the  dimensionless  notation 


.  Z  £  U  4  V  V 

W  -  -  .  5  =  —  ,  4  =  —  ,  T  =  T  (  s) 


a2  («) 


U, 


u, 


(Bl) 


(B2) 


where  t(£)  is  the  relative  temperature  distribution,  expressed  as  a  function 
of  the  relative  velocity,  Kj  becomes 


K,  * 

where 


sec  al(ec).  <B3) 


- (B4a) 

< e  •  ec> 


t<£)  c'tfeHS  ,  (B4b) 

<e  -  se'4 


and  we  have  assumed  U  =  (  £  =  1 )  at  z  =  5  ( £  =  1 ) . 

Integrating  Eq.  (B4b)  once  by  part%  along  a  path  indented  over  the 

singularity  at  4  =  4  .  and  then  separating  out  the  singular  portion  (I  ), 

C  8 

we  obtain 


I  *  -T(0)t'(0)6;1  -  4'(U(i  -  $ J"1  +  I  +  [  — !T-;  -c-fo 


<e-ec) 


(B5) 
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where 


(B6) 


and  the  subscript  c  implies  evaluation  at  the  singular  point  £  =  £c . 
Integrating  I  over  the  singular  point  yields  directly 

s  — —  -i 


(B7) 


Substituting  this  result  in  Eq.  (Bl)  yields  for  the  imaginary  part  of  Kj 


ImKj  =  -  sec2a  (rC)^ 


<B8a) 


irMj  sec  06 


1 

[ufl 

A 

~  T 

ds 

U*  <«>_ 

(B8b) 


a  sb 


e  ~  1  2  2„ 

to  a  sec  a 

c 


T<«c) 

U,3C«c) 


A  [~U'  <»>' 

ds  Lt  (*). 


(B8c) 


s=z 


We  consider,  as  a  particular  example,  the  turbulent  velocity  profile 


;  =  g 


(B9) 


together  with  the  quadratic  (in  velocity)  temperature  profile 

T 

9  ^  =  _f3Lr_M  w 

1 


t  ( t )  *  a  v  b£  +  c  £*,  a  =  —  ,  c  = 

T, 


M?  ,  b  -  1  -  a  -  c, 


(BIO) 


1 


for  which 


(r;’)'  =  7  (6  a  £5  +  7b  £6  +  8c  {). 


(Bli) 
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We  not*  that  Eq.  (B9)  is  not  accurate  right  down  to  l  »  0,  where  it  gives 
t '  (0)  ■  0,  hut  the  error  in  equating  the  first  term  in  Eq.  (B5)  to  sero  for 
large  free  stream  Reynolds  numbers  <say  greater  than  10  )  is  negligible. 
Substituting  Eqs.  (B9)  through  (Bli)  in  Eq.  (B5)  and  carrying  out  the 
integrations  yields 

1  -  ?[al5(€c)  +  bl6(fic)  +  cl?(«c)  -  (1  -  ec)“4J  .  (B12) 


where 


n-1 

I  (t  )  .  (n  ♦  t)i  S  ♦  i 

[_m=0  n  -  m 


(B13) 


Substituting  Eq.  (B12)  in  Eq.  (B3)  and  a,  b,  and  c  from  Eq.  (BIO)  then  yields 
K,  •  -  1  +  7MjZ  M»c2  «j[l6<ec>  -  o  -  c.)-1] 


+Mm*  [w-wj 


The  quantity  actually  required  in  the  stability  calculations  is 
mflt,  =  j^l - Kj  *  M2  cos2a-(l  -  , 


(B14) 


(B15) 


which  we  write  in  the  form 


mf  K, 


A<ec)  +  B(ec)  +(Y  * 1 


2  ,Mr  C'«c» 


(B16) 


with 


A(4e)  •  7(1  -{/  [x5(tc)-l6(tc)] 


(Bi  7) 
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Bit.)  »  7(1  -  tc)  [(1  -  !c)  I6(tc)  -  lj  (B18) 

C(«c)  *  7(1  -Sc)2[l6(ic)-I7(Sc)]  •  <*!9> 

The  real  and  imaginary  parts  of  A,  B,  and  C  are  plotted  in  Figs.  7a,  b,  c, 
respectively.  The  real  and  imaginary  parts  of  m^ are  plotted  in 
Figs,  la,  b  for  a  =  0,  M  =0.1,  1.2,  1.5,  and  2.0,  and  an  insulated 
boundary- -i.  e.  ,  one  for  which 


T 


o 


T 


1 


(B20) 


The  latter  results  also  may  be  used,  with  only  a  small  error,  for  a  ^  0  if 
the  nominal  value  of  Mj  on  the  curves  is  replaced  by  cos  a. 


i  « 


Figure  7c,  C  (£),  at 
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|  APPENDIX  C 

MAXIMUM  VALUE  OP  Im  v 

We  wish  to  maximise  the  quantity 


-  Im  v  *  ( \JT/Z)  A-B  •  A  +  B  -  2  0  j 

(Cl) 

with  respect  to  £,  A  and  B  being  given  by 

H 

r  3  T  32  2  3  I 

'  (1  -  3  £0  4  0A)  a  (1  -  3  £0  +  eV  -  (£  +  0  )  J  J 

(C2) 

Differentiating  Eq.  (Cl)  with  respect  to  £  and  equating  the  result  to 
sero,  we  obtain 

Q(A  +  B)8  +  ABj  (A  +  B  -  29)  a  0. 

(C3) 

It  is  evident  that  the  maximum  value  of  -Im  v  corresponds  to  the  vanishing 
of  the  first  factor  in  Eq.  (C3);  dividing  this  through  by  ABO  and  combining 
the  resulting  equation  with  the  values  of  A  B  and  A3  +  B3,  obtained  from 

Eq.  (Cl),  we  have 

a*1  +  b"1  +  e"1  =  0, 

(C4a) 

and 

AB  -  (£  +02)  =  0, 

A3  +  B3  -  2  (1  -  3£6  +03)  =  0. 

(C4b) 

(C4c) 

Eliminating  A  and  B  from  Eqa.  (C4a,  b,  c)  yields  the  cubic  equation 

£3  -  904£  +  203  «  0  (C5) 

for  the  required  value  of  £,  say  £m  .  The  corresponding  value  of  -Im  v 
may  be  placed  in  the  form 

(Jmv)max  *  31/2  ’  2’1/3Mn>.  n  •  31/20  .  (C6) 

where  h,  normalised  to  have  the  value  unity  at  t|  ■  0,  is  given  by 

h  «  •  1 1  +  1*  |  •  [u  -  ’!*)  (  t  +  -J-  1/2  .  (C7) 
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and  r\  is  introduced  in  place  of  0  bacauaa  of  tha  manner  in  which  it  enter* 

t  (^i.). 


£q,  (C5)  haa  only  one  real  root,  opposite  in  sign  to  0,  if  j  q  j  ■  |  3l^*Q  |  <  i, 
while  if  j  n  j  >1  there  are  three  real  roots;  solving  for  these  roots  yields 


If  |n  >  1  that  value  of  {  yielding  the  maximum  value  of  h  is  to  be  chosen. 


The  limiting  values  of  $  for  large  i  r\  I  are 


£  s*n|n  -  0-3/2  -1  ,.-3/2-1 

»m  —  T1  |  *1  m  *  *1  ,  6  •  J  'l. 


The  corresponding,  asymptotic  forms  of  h  are 


h—  2 


,4/3  .-1/4  „  1/2 


n  »i 
■T)  » 1, 


(CIO) 


(Cl  la) 
(Cl  lb) 


We  note  that  these  approximations  actually  are  quite  adequate  for 


r\  >1. 
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